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Finite temperature Casimir effect for charged massless scalars in a magnetic field 
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The zeta function regularization technique is used to study the finite temperature Casimir effect 
for a charged and massless scalar field confined between parallel plates and satisfying Dirichlet 
boundary conditions at the plates. A magnetic field perpendicular to the plates is included. Three 
equivalent expressions for the zeta function are obtained, which are exact to all orders in the magnetic 
field strength, temperature and plate distance. These expressions of the zeta function are used to 
calculate the Helmholtz free energy of the scalar field and the pressure on the plates, in the case 
of high temperature, small plate distance and strong magnetic field. In all cases, simple analytic 
expressions are obtained for the free energy and pressure which are accurate and valid for practically 
all values of temperature, plate distance and magnetic field. 
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I. INTRODUCTION 



The Casimir effect is a quantum phenomenon where an attractive or repulsive force is observed between electrically 

neutral conducting plates in vacuum, and can be regarded as a quantitative proof of the quantum fluctuations of the 

electromagnetic field. Casimir first discovered the effect, by calculating the attractive electromagnetic force between 

two parallel conducting plates [1|. The repulsive Casimir effect was discovered by Boyer some time later, when he 

showed that if the electromagnetic field is confined inside a perfectly conducting sphere, the wall of the sphere is 

subject to a repulsive force [2j. The first experimental evidence of the Casimir force was obtained more than fifty 

years ago by Sparnaay [3| and, since then, many greatly improved experimental observations have been reported. For 

I ■ a comprehensive review of these experiments, see the review article by Bordag et al. 4J. 

^. | Since Casimir forces have many applications -from nanotubes and nanotechnology |5|-l8| , to branes and compactificd 

£"^ i extra dimensions J9H28I] . to string theory |29h32j- a large effort has gone into studying the Casimir effect and its 

t-H ' generalization to quantum fields other than the electromagnetic field: fcrmions 33], bosons and other scalar fields 

have all been investigated extensively J4j]. 
^^ ' It is well known that Casimir forces are very sensitive to the boundary conditions of the involved quantum fields on 

the plates. In the case of scalar fields, the most used boundary conditions are Dirichlet and Neuman, in the case of 
fcrmion fields or fields with spin in general [34J bag boundary conditions are used. In this work we will use Dirichlet 
C*i | boundary conditions for a scalar field confined between two parallel plates. 

Scalar fields, with or without charge or mass, appear in many different areas of physics: the Higgs field is responsible 
for spontaneous symmetry breaking in the Standard Model, dilaton fields break the conformal symmetry of the string 
in superstring theory [351 136 |. inflatons are used to solve the problem of a non- vanishing cosmological constant by 



causing the accelerated expansion of the universe [37H39J . In condensed matter physics, scalar fields are important 
to describe spontaneous breaking of discrete symmetries: the Ginzburg-Landau scalar field is associated with type II 
superconductors, scalar fields are used to explain Landau diamagnetism [4CI l4l|. etc. 

The Casimir effect for charged scalar fields in a magnetic field has been studied in vacuum [42| and at finite 
temperature [43| using the Schwinger proper time method to calculate the effective action, but these authors are only 
able to obtain the free energy as an infinite sum of modified Bcssel functions. In this paper we use a different method, 
the zeta function technique, to study the Casimir effect for massless scalar fields at finite temperature and in the 
presence of a magnetic field. This method allows us to obtain simple analytic forms for the free energy and Casimir 
pressure, valid for practically all values of the parameters involved. A similar investigation of the Casimir effect for 
massive scalar fields at finite temperature and in the presence of a magnetic field will be presented elsewhere. 

Since Casimir effect calculations must carefully address the issue of regularizing the vacuum energy, it is best to 
use the most effective regularization techniques. Many regularization techniques are available nowadays, and two of 
them have been applied freq uently and successfully to the Casimir effect: the cutoff method often used in various 
piston configurations [4J, |45| and the zeta function technique [4q448| . As we stated above, the choice for this paper 
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is the zeta function technique, a powerful regularization technique used also in the computation of effective actions 
49, 50]. We apply this regularization to obtain the free energy and Casimir pressure due to a scalar held confined 
between two parallel plates, at a distance a from each other. We assume Dirichlet boundary conditions on the plates 
and take our system to be in thermal equilibrium with a heat reservoir at finite temperature T, using the imaginary 
time formalism of hnitc temperature field theory, which is suitable for a system in thermal equilibrium. A uniform 
magnetic field B is present in the region between the plates and is perpendicular to the plates. 

In Section[nl we obtain three equivalent expressions of the zeta function for this system, exact to all orders in eB, T 
and a, where e is the scalar field charge. We also obtain simple analytic expressions for the zeta function in the case of 
high temperature (T ^ a -1 , \fe~B), small plate distance (aT 1 ~^> T, y/eB), and strong magnetic field {-/eB ~^> a~ x ^T) 
In Section UITl we use the zeta function obtained in the previous section, to calculate the Helmholtz free energy of the 
scalar field and the pressure on the plates and obtain simple analytic expressions for these quantities in the case of 
high temperature, small plate distance, and strong magnetic field. A discussion of our results is presented in Section 

El 

II. ZETA FUNCTION EVALUATION 

Using the imaginary time formalism of finite temperature field theory, we write the partition function Z for a 
bosonic system in thermal equilibrium at finite temperature T 

Z = N I Dc/)* D(j> exp [ I dr I d 3 xC J , (I) 

./Periodic \J J J 

where C is the Lagrangian density for the bosonic system, N is a constant and 'periodic' means that this functional 
integral is evaluated over field configurations satisfying 

c/)(x,y,z,0) =<j>(x,y,z,P), (2) 

where j3 = 1/T is the periodic length in the Euclidean time axis. In addition to the hnite temperature boundary 
conditions given by @, we impose Dirichlet boundary conditions for scalar bosons between two square plates. In 
3-dimensional space with two large parallel plates perpendicular to the z-axis and located at z = and z = a, the 
Dirichlet boundary conditions constrain the scalar field to vanish at the plates, 

(f>(x,V,0,T~) = 4>(x,y,a,T) =0. (3) 

In the slab region there is also a uniform magnetic field pointing in the z direction, B — (0, 0, B). The scalar field has 
charge e and thus will interact with the magnetic held. 

The scalar held Helmholtz free energy F and partition function Z are related by 

F = -/T 1 logZ. (4) 

A straightforward evaluation of the functional integral ([1]) yields 

logZ = -logdet(- J D E |^ a ), (5) 

where the symbol T a indicates the set of functions which satisfy boundary conditions ([2]) and ([3]) , and the operator 
-De is defined as: 

D^ = d 2 + d 2 -(f-eA) 2 ± , (6) 

where A is the electromagnetic vector potential, the subscript E indicates Euclidean time, and we used the notation 

P± = (Px,Py,0). 

The zeta function technique allows us to use the eigenvalues of De to evaluate logZ. The Dirichlet boundary 
conditions ([3]) are satisfied only if the allowed values for the momentum in the z-direction are 

7T 

Pz = -n, (7) 

a 

where n € {0, 1, 2, 3, ...}, and therefore the eigenvalues of —d 2 . — d 2 whose eigenfunctions satisfy © and §S§ are: 

7T 2 9 47T 2 9 



where n G {0,1,2,3,...} and m G {0, ±1, ±2, ±3, ...}. The spectrum of the operator (p — eA)\ is well known from 
one-particle quantum mechanics, and its eigenvalues are the Landau levels 



2eB I 



1 
2) ' 

with I G {0, 1, 2, 3, ...}. Using the eigenvalues ([5]) and ©, we construct the zeta function £ (s), which is given by 
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where L 2 is the area of the plates. Once we put ( in a suitable closed form, using the zeta function technique we will 
obtain immediately the partition function 



and then the free energy using (HJ) 
With the help of the following identity 



logZ = C'(0), 



r( s ) J 

where T(s) is the Euler gamma function, we rewrite £(s) as 
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where we also used 
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2 sinh z 



(11) 
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(13) 
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It is not possible to evaluate (1141) in closed form for arbitrary values of B, a and /3, but it is possible to find simple 
expressions for £(s) when one or some of B, a and T are small or large. From these simple expressions of the zeta 
function, the free energy will be obtained immediately. 

First we evaluate £(s) in the high temperature limit. To do so, we apply the Poisson resummation formula [5l| to 
the n-sum in (Tl4l) and obtain 



CO) = a (b(s) + Cs,a(«) + (B,T(s) + CB.a.T(s) 



where 
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(19) 
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After changing the integration variable from t to tnaf3/2irm in (|20l) , we obtain 
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(21) 



When 2aT 3> 1, only the term with n = m = 1 contributes significantly to the double sum so, using the saddle point 
method, we evaluate the integral for eB <C 4ir 2 T 2 and obtain 



Cb,o,t(s) = 



L 2 eB /ojSV e" 4 ™^ 



27rar(s) V2^7 sinh(^ 



-eBa0\ 



Next we evaluate (|18p for eB <C 47r 2 T 2 . In this case, we can set 

— « 1 - -(eBt) 2 

sinh eBt 6 V y 

and, after substituting (|23|) into (fl"8|) , we integrate to find 



Cb,t(s) 



L 2 ( /3 



T(s) \2tt 



2tt 3 / 2 
/3 3 
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(22) 



(23) 



(24) 



where Cr is the Riemann zeta function of number theory. For calculating the free energy, we only need to know £(s) 
for s — > 0. For small s we have 



z 2s ( R (2s - 3) 



r(s - |) _ y¥ 
r( s ) 90 



s + G(s 2 ), 



and 



^CH(2 S + l)^^ = ^ + v^(7 B + ln0.s + O( S 2 ), 



where 7e = 0.5772 is the Euler Mascheroni constant. Substituting ([25]) and (|26|l into (J24J) , we obtain 

r 2 e 2 B 2 /?.l , /3 
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(25) 



(26) 



(27) 



valid for ei? <§C 47r 2 T 2 and small s. 

In the high temperature limit, both eBa 2 -C 1 and eBa 2 3> 1 are possible, and therefore we need to evaluate (fig)) 
for both scenarios. When eBa 2 Clwe use (|23|) into (fi~9|) . integrate, and find 
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which, for small s, becomes 



L 2 a 



47rV 2 r(s) 



V(| - s )Ca(3 - 2s) - ^-^r(-i - a )C fl (-l - 2a) 



() 



(28) 



Cs,a(s) 



where 0?(3) = 1.2021. When eBa 2 > 1 we use 
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into (|T9|) . change the integration variable from £ to -9%tj and find 



Cs, a (s) 
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(29) 
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Only the term with n — 1 contributes significantly to the sum when eBa 2 3> 1 and, using the saddle point method, 
we evaluate the integral and find 
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Finally we calculate Cb(s), the only piece of the zeta function that can be evaluated exactly and, after integrating, 
we find 



r2( p p>\3/2-s 

Cs(s) = 4 U»r(«) (1 - 2V " s)r(s - ^ Cr{s - ^ 



(33) 



which, for small s, becomes 



T2( f>\3/2 

Us) = { J (72- l)Cfl(-i)*, 



(34) 



where &*(-*) = -0.2079. 

By adding fl22]), P?]). (|2"9"|) and ([51],. we find £(s) in the high temperature and very weak field limit, 2T > a -1 , 
IT > y/eB/ir, and eB < a" 2 
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where we took the small s limit. By adding (1221) . (|27|l . (|32|) and (l34|) , we find £(s) in the high temperature and very 
large plate distance limit, 2T > a~ l , 2T > \feBj-K, and e/3 > a -2 



C(«) - i 2 



7r 2 a (eBf/ 2 a r r SJ _ . . eB e r^ a IP eB 2%/77? „ e 2 B 2 (3aA , /3 , 

45^ + HH^" ^-*) + ^^^) + 2^ ~ "isf^ + ^ + ^\ 



s, (36) 



where we also took the small s limit. 

Next we evaluate £(s) in the limit of small plate distance and apply the Poisson resummation formula to the m-sum 
in (TT4"|) to obtain 



C(*) = ~ |C/;iM --,, i,-,,i>i -.- ku. 7 1*)+^. .,.,(«) 



where Cb( s ) is the same as in (JTTJ), and 

Cs,a(s) = 



r2 °° /-oo 



47r 3 / 2 r( s ) ^ J 
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(38) 
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It is evident from (J3TJ) - (gDJ) that £(s), in the limits 2aT < 1 and e_B < 7r 2 a~ 2 , is obtained from fll6) - ([20]) by 
replacing a with /3/2 and (3 with 2a. For eB ( | J <1 and small s, we find 



C( S ) = L 2 
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for eB ( £ ) 3> 1 and small s we find 

C( S )=L 2 
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Last, we evaluate £(s) in the strong magnetic field limit, eB ^ ( f ) anc ^ e -^ ^ a 2 - Under these conditions, 
after applying the Poisson resummation formula to both the n- and m-sum in (|14[) , we find 



((s)=al3[Cw(s) + as)} 



where 
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is the zeta function of the one-loop vacuum effective lagrangian for massless scalar QED first calculated by Weisskopf, 
and 
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The integral in (PHI) can be evaluated exactly, and we find 

L 2 (eB) 2 ~ s 
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which, for small s, becomes 
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where we used the interesting numerical fact [55 
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We evaluate £(s) by using (|30l) . which is valid in the strong magnetic field limit, then change integration variable 
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where the term with m = n = is excluded and only terms with n = 0, ±1 and m = 0, ±1 contribute significantly to 
the double sum. We integrate using the saddle point method and, for small s, obtain 
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Adding (l4"Tl) to (I5TH) we find the zeta function in the limit of strong magnetic field, eB ^$> a 2 and eB ^> ( § ) > 
small s 



and 
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III. FREE ENERGY AND CASIMIR PRESSURE 

The derivative of the zeta function is obtained easily by taking advantage of the useful fact that, for a well behaved 
G(s), the derivative of G(s)/T(s) at s = is simply G(0) [53[ and therefore, using ([T2")) and ([Till , we find the free 
energy 
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Using our results for the zeta function ([16]) , (|37[) . and (|43|) . we are able to obtain three other expressions of the free 
energy, all equivalent to (J52)) . 
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best suited for high temperature expansion (2Ta ^> 1 and 2T 3> y/eB/n) 
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best suited for small plate distance expansion (2Ta -C 1 and a 1 ^$> s/eB/n), and 
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best suited for strong magnetic field expansion. The last equation has been obtained by other authors [43J , who used 
(1551) to write the free energy as an infinite sum of modified Bessel functions. 

It is not possible to evaluate (|52j) - (|55|) in closed form for arbitrary values of B, a and /3 but, using our results 
from Sec. HJ we found simple analytic expressions for the free energy when one or some of those three quantities are 
small or large. To calculate the free energy in the high temperature limit, 2T » a -1 and 2T » y/eB/n, we use ((35)) 
and (1361) to find 
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valid for eB 3> a 2 . Notice that in ([55)) and ([571) the dominant term is the Stefan-Boltzmann term — ^rVT 4 , where 
V = L 2 a is the volume of the slab. Terms with a linear dependence on the plate distance, such as this one, are 
proportional to the volume of the slab and represent a uniform energy density. If the same magnetic field is present 
outside the slab and the medium outside the slab is also at temperature T, such terms do not contribute to the 
Casimir pressure. If there is vacuum outside the slab, i.e. no magnetic field and zero temperature, uniform energy 
density terms contribute a constant pressure which is very easily calculated. In this paper we assume that the same 
magnetic field is present inside and outside the slab, and that the medium outside the slab is at the same temperature 
as the one inside the slab, so we neglect contributions to the Casimir pressure from uniform energy density terms. 
The pressure P on the plates is given by 
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and therefore, for 2T 3> a 1 ^> y/eB/ir, we find 
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and 
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for 2T S> s/eB/n S> a x . Since the third term in (|59|) is negligible when compared to the other ones in ([59| and (|60[) , 
we write the high temperature Casimir pressure as 
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for2T> ^B/vr>a- 1 . 

Next we obtain the free energy in the small plate distance limit, using (f4"Tj) and (|12"1) . to find 
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for a 4 3> \/eB/TT 3> 2T. The dominant term here is ~^npr, which is the familiar vacuum Casimir energy for a 
complex scalar field, and for the photon field QJ. The Casimir pressure for small plate distance is 
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in the case of very weak magnetic field (a 1 3> 2T ^> y/eB/ir), and it is identical in the case of very low temperature 
(a -1 3> y/eB/n 3> 2T). Since the third term in (|65[) is much smaller than the other ones, we can neglect it and write 
the pressure in the small plate distance limit as 



eB 



-■K/3/a 



' 240a 4 + 2^ S inh(^ff£) ~ 48tt 2 ( n 2^ 



1) 



Finally, for strong magnetic field the free energy is found using (f5~Tj) 



F = -L 2 a 



2d2 / i 

-Ar lneS-ln3-- 

96tt 2 V 2 



(eB) 



5/4 



,-2aVeS 



2tt 3 / 2 



2a 3 / 2 



£3/2 



-2VeBy/a 2 +0 2 /4^ 



(a 2 + /3 2 /4) 



3/4 



(66) 



(67) 



where the dominant term is the one-loop vacuum effective potential for massless scalar QED [52J , and is proportional 
to the volume of the slab, as expected. The effective potential is a uniform energy density term and therefore, under 

our assumptions, does not contribute to the Casimir pressure. The pressure, for eB 3> ( £ 
given by 



and eB ^> a 2 , is 



P = 



(eB) 5 / 4 



-2aVeB 



'eB 



a 



2tt 3 / 2 ^ 
and, neglecting the smaller terms, we obtain 

(eB) 7 / 4 



— 1_\ (eB) 5 / 4 e ^^V° 2 +^/4 / 2q 2 

eB+ ^) + 2vr 3 / 2 (a 2 + ^2/4)3/4 ^ " V^+ /?2/4 2a 2 + /3 2 /4 / 



r 3/2 



-2aVeB 



2y/E 



a 2 e -2V^By/a2+{3y4 S 

(a 2 +/? 2 /4) 5/4 



(68) 



(69) 



IV. DISCUSSION AND CONCLUSIONS 

In this paper we used the zeta function regularization technique to study the finite temperature Casimir effect of 
a massless charged scalar field confined between parallel plates and in the presence of a magnetic field perpendicular 
to the plates. We have obtained three expressions for the zeta function (Tip]) . (j3"T)l . and (H51) . which are exact to 
all orders in the magnetic field strength B, plate distance a and inverse temperature /?, and used them to derive 
expressions for the Hclmholtz free energy and for the Casimir pressure on the plates in the case of high temperature 
(AT 2 > a -2 , eB/ir 2 ), small plate distance (a~ 2 > AT 2 , eB/ir 2 ) and strong magnetic field (eB > a~ 2 , AT 2 ). 

We have been able to numerically evaluate the free energy with very high precision, using the three exact expressions 
(f53|) . ([54]) . and (f55|) . and we compared the values of the free energy obtained from our simple analytic expressions to 
the exact numerical values. In the high temperature case we found that, for 2aT = A, Eq. (156j) is within 0.7 percent 
of the exact value of the free energy in the range < eBa 2 < 1, while Eq. (|57|) is within 0.7 percent of the exact 
value of the free energy in the range 1 < eBa 2 < oo. For 2aT — 10, Eq. (I55|) is within 0.05 percent of the exact 
value of the free energy in the range < eBa 2 < 1, while Eq. (|57|) is within 0.05 percent of the exact value of the 
free energy in the range 1 < eBa 2 < oo, showing a very rapid decrease of the small discrepancy between Eqs. (1561) . 
(|57[) and the exact values of the free energy. We summarize this finding by writing the free energy per unit area in 
the high temperature limit as 

f ^ a (^-lRfl(-JQ(e-B) 3/2 ° eB e -«"//i CfiO) , e'aV , e 2 B 2 g s , , s f nr n<„n„2<,. 

F \ ~IEm 2^3 2^ sinh (£|^y 8^3^ + 96tt/3 + TSvr^ W E + in 17 ) t0r U - etia - 1 ' 

L 2 ~ I _jSa _ (V2-lKn(-^(eBf^a _ eB e -^'"p _ eB_ 2 ^Ba £&£(„„ , K, A) fo r 1 < eB a 2 < OO 

[ 45^ 277/3 277/3 s inh(iM^Zy 2tt/3 e + 48t J 17S + m 47r I IOr 1 ^ e£SO <00. 

Equation (|70[) is a simple analytic expression of -F in the high temperature limit, valid for all values of the magnetic 
field B and the plate distance a, and with a discrepancy of no more than 0.7 percent from the exact value of F for 
2aT > 4. A similarly accurate expression of the Casimir pressure P, valid for 2aT > A and all values of a and B, is 
obtained immediately from ((70)1 , since P = —j^^K;- 

In the small plate distance case we found that, for 2aT = |, Eq. (1551) is within 0.7 percent of the exact value of the 

free energy in the range < eB I ^ j < 1, while Eq. (|64|) is within 0.7 percent of the exact value of the free energy 
in the range 1 < eB I ^ ) < oo. For 2aT — j^, Eq. ([55)1 is within 0.05 percent of the exact value of the free energy 
in the range < eB I g ) < 1, while Eq. (|M|) is within 0.05 percent of the exact value of the free energy in the range 

1 < eB I ^ ) < oo, showing again a very rapid decrease of the small discrepancy between our analytical expressions 

and the exact values of the free energy. We summarize the small plate distance limit by writing the free energy per 
unit area as 

F I 720^ 4tt 2tt/3 ci nh( "f f ) ^W + 384tt + ^8^ W B + m 2n > IOT U - CD \ 2 ) - L ' 

[ 720?- 4^ 2^ sinh( £BM7~ 27p e + ^8^-^E + m^) tor 1 < e±J (^ ^ J < OO, 

(71) 
a simple analytic expression of F, valid for all values of B and T, and with a discrepancy of no more than 0.7 percent 
from the exact value of F for 2aT < j . The pressure in the case of small plate distance is obtained immediately from 
(f7i"j) for 2aT < i and all values of B and T. Notice that, if we set T = in (|7"Tj) . we obtain the Casimir energy Eq 
for a massless and charged scalar field in a magnetic field, 

Ec n 2 (V2-lK R (-i)(eBf/ 2 e 2 B 2 a ( a 

lJ = -t2^ T, + l8^ (7£ + ln ^ } ' (72) 

where we see that the magnetic field, as it grows, inhibits the Casimir energy of the scalar field j42[. Our result, a 
simple analytic expression for Ec , is more explicit than that of |42J | where the magnetic field correction to the Casimir 
energy is presented as an infinite sum of integrals. 

In the case of strong magnetic field, the free energy shown in Eq. (|67|) is valid for all values of a and T, and so is 
the pressure shown in Eq. (|69p . If we set T = in (I67|) . we can neglect the effective potential which is a uniform 
energy density term, and obtain Ec in the strong magnetic field case 

Ec 3 (eBfl^e- 2 ^ 

L 2 47r3/ 2 JE [ ' 
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which agrees with [42[ on the dependence of Ec from a and B, but is in disagreement for the overall sign, since we 
obtain a negative value for Eq, not a positive one. We also obtain the numerical constant present in Ec, while the 
authors of Ref. 421 did not. 
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